The pinch technique at two-loops: 
The case of mass-less Yang-Mills theories 

Joannis Papavassiliou 

Departamento de Fisica Teorica, Univ. Valencia 
E-46100 Burjassot (Valencia), Spain 



ABSTRACT 

The generalization of the pinch technique beyond one loop is presented. It is shown that the crucial physical principles 
of gauge-invariance, unitarity, and gauge-fixing-parameter independence are instrumental for accomplishing this task, 
• and is explained how the aforementioned requirements single out at two loops exactly the same algorithm which has 
' been used to define the pinch technique at one loop, without any additional assumptions. The two-loop construction of 
0^ . the pinch technique gluon self-energy, and quark-gluon vertex are carried out in detail for the case of mass-less Yang- 
I Mills theories, such as perturbative QCD. We present two different but complementary derivations. First we carry 
O . out the construction by directly rearranging two-loop diagrams. The analysis reveals that, quite interestingly, the 
2^ well-known one-loop correspondence between the pinch technique and the background field method in the Feynman 
gauge persists also at two-loops. Since we use dimensional regularization the entire construction does not depend 
on the value of the space-time dimension d. The renormalization (when d = A) is discussed in detail, and is shown 
to respect the aforementioned correspondence. Second, we present an absorptive derivation, exploiting the unitarity 
of the S'-matrix and the underlying BRS symmetry; at this stage we deal only with tree-level and one-loop physical 
amplitudes. The gauge-invariant sub-amplitudes defined by means of this absorptive construction correspond precisely 
^Jq I to the imaginary parts of the n-point functions defined in the full two-loop derivation, thus furnishing a highly non- 
. trivial self-consistency check for the entire method. Various future applications are briefly discussed. 
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I. INTRODUCTION 



The pinch technique (PT) was introduced by Cornwall [||, as the first step in the development of a self-consistent 
truncation scheme for the Schwingcr-Dyson equations of QCD It was further developed in and was generalized 
to theories with spontaneous symmetry breaking , and subsequently to the full Standard Model . Various phe- 
nomenological applications were presented M, and the one- loop correspondence between the PT and the background 
field method (BFM) [p|~p^ was established . The important role of unitarity and analyticity was appreciated and 
exploited in a series of papers |[l^-|l5t , which allowed the generalization of the Breit-Wigner formalism for unstable 
particles to the non-abelian sector of the Standard Model. Parallel developments took place in the field of finite 
temperature QCD the concept of the QCD effective charge Q was examined in detail [0, and subsequently 
applied to the field of renormalon calculus , and the generalized PT was presented . In addition, a formalism for 
studying resonant CP violation ||2^ has been developed , and various issues related to the high-energy behaviour 
of the W- fusion process have been resolved Q . 

The PT reorganises systematically a given physical amplitude into sub-amplitudes, which have the same kinematic 
properties as conventional n-point functions, (propagators, vertices, boxes) j2^, and are in addition endowed with 
important physical properties. The essential ingredient one uses in this rearrangement is the full exploitation of the 
elementary Ward identities (WI) of the theory, in order to enforce crucial cancellations. The various types of physical 
problems for which the PT can serve as a useful tool have been discussed frequently in the literature, most recently 
in (H] and [l^. 

So far the PT program has been carried out only at one-loop order, and its generalization to higher orders has 
been a long-standing question . The general methodology of how the extension of the PT proceeds at two-loops 
has been presented in a recent brief communication |p6| . There it was explained how the one-loop formalism, when 
properly interpreted and suitably adopted to the two- loop context, leads naturally to the PT extension. In this paper 
we will address the technical aspects of this procedure in detail, and will discuss extensively the plethora of physical 
issues involved. 

We will present two independent constructions, which arc however inextricably connected. The first one deals 
directly with the two-loop ^-matrix element for the scattering process qq qq, exactly as was the case in the early 
one-loop PT applications. The first crucial ingredient is that of gauge-fixing parameter (GFP) independence. 5*- 
matrix elements are guaranteed to be GFP independent; in fact, as happens in the one-loop case, the cancellation 
of all GFP dependent terms proceeds by exploiting tree-level WI only, without having to carry out sub-integrations. 
This property is important, because it preserves the diagrammatic representation of the S-matrix, as well as the 
kinematic identity of the sub-amplitudes appearing in it (propagators, vertices,boxes). One can therefore, without 
loss of generality, choose a convenient gauge, such as the renormalizable Feynman gauge (RFG), provided that one 
considers the entire set of two-loop Feynman graphs contributing to the given S'-matrix element 

In the one-loop case the next step would be to split the bare three-gluon vertex appearing inside the one-loop 
quark-gluon vertex into a pinching and non-pinching contribution. This splitting is very special because it guarantees 
that the resulting effective Green's functions satisfy (at one-loop) naive, QED-like WI Q, instead of the usual Slavnov- 
Taylor identities As a consequence the effective one-loop gluon self-energy captures the running of the QCD 

coupling, exactly as happens with the photon vacuum polarization in QED. In the two- loop case the construction is 
lengthier but the crucial operation is precisely the same. One carries out the aforementioned splitting to all vertices 
whose one of the incoming momenta is the physical momentum transfer (or center-of-mass energy) of the process. As 
we will see in detail this splitting is sufficient to give rise to two-loop PT effective Green's functions which have the 
exact same properties as their one-loop counterparts. 

The reason why all other three-gluon vertices inside the loops should remain unchanged (no splitting) can be 
understood by resorting to the special unitarity properties that the PT sub-amplitude must satisfy, a point which 
brings us to the second derivation. For this derivation we employ the unitarity and analyticity properties of the 
S'-matrix, very much in the spirit presented in the more recent one-loop PT literature There, the precise 

diagrammatic correspondence between the one-loop (forward) process qq — > qq and the tree-level process qq — > gg 
provided very useful, stringent constraints on the entire construction, rendering the method all the more powerful. 
These constraints are automatically encoded in the two-loop PT construction presented here; in fact they are even 
more constraining than in the previous order, for reasons that we will describe qualitatively now, and in great detail 
in the main body of the paper. The imaginary parts of the two-loop PT Green's functions are related by the optical 
theorem to precisely identifiable and very special parts of two different on-shell processes, the one-loop process qq gg 
and the tree-level process qq ggg. In particular, the two-particle Cutkosky cuts of the two-loop PT self-energy 
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are related to the PT rearranged s-channel part of the one-loop qq gg, while, at the same time, the three-particle 
Cutkosky cuts of the same quantity are related to the PT rearranged s-channel part of the tree-level qq ggg; the 
latter is the exact analogue of the PT rearranged s-channel part of the tree-level qq gg, already studied in the 
literature cited above. 

The paper is organized as follows: In section II we present a brief overview of the one-loop PT algorithm, and 
discuss the most characteristic properties of the one-loop PT effective Green's functions. Here we follow the original 
PT formulation |^,||, and postpone the overview of the one- loop absorptive PT construction [^,0 until section 
V. In section III we present the full two-loop construction. This section contains three sub-sections, where we deal 
separately with the one-particle reducible graphs, the two-loop quark-gluon vertex, and the two-loop gluon self-energy. 
In the second sub-section we use at intermediate steps various results for the one-loop three-gluon vertex which are 
derived in section VI, and are more naturally integrated in the analysis presented there. Section III contains the main 
results of this paper; a major highlight is the proof that the correspondence between the PT and the Background Field 
Method Feynman gauge (BFMFG) persists at two-loops. In section IV we carry out the renormalization program in 
detail. We show that the two-loop PT Green's function can be renormalized by judiciously rearranging the existing 
counterterms, and that this procedure respects the correspondence between PT and BFMFG established in the 
previous section. In section V we present the general formalism and methodology of the absorptive PT construction; 
this second derivation is completely independent of the first, but at the same time is deeply connected to it. In this 
section we first derive various formulas which will be used in the next section, and present a thorough overview of 
the one-loop case. In section VI we present the two-loop absorptive derivation. There are two sub-sections: the first 
one contains a detailed discussion of the one-loop process qq — *■ gg, and its role in enforcing the unitarity of the 
individual two-loop PT Green's functions; the second sub-section contains the study of the PT rearranged tree-level 
process qq — + ggg, which demonstrates the crucial function of this process in realizing the aforementioned unitarity 
properties. Finally, in section VII we present our conclusions, and discuss possible connections with other work as 
well as future applications. 



II. OVERVIEW OF THE ONE-LOOP CASE 



In this section we briefly review the one-loop PT construction and establish some useful notation. 
The fundamental tree-level three-gluon vertex Tal_iy{q,pi,p2) is given by the following manifestly Bose-symmetric 
expression (all momenta are incoming, i.e. q + pi + P2 = 0) 

raji/(9,Pi,P2) = {q~pi)ugat, + {pi -P2)ag^,u + {p2 - q)^,gau ■ (2.1) 

^alviq,Pi,P2) may be split into two parts 

r^°Uq,Pl,P2) = rP^,Aq,pi,P2) + rP^^Aq,pi,P2) , (2.2) 

with 

ri?i^i.('7,Pl,P2) = {Pl -P2)agf,^ + "^q^gaf, - 2(7p5ai/ , 

^plf,,yil^Pl^P2) ^ P2ugaf, - Plfj^ga,. ■ (2.3) 

The vertex ^^p^^^iq, Pi, P2) is Bose-symmetric only with respect to the fi and ly legs. The first term in P^^^^^ is a 
convective vertex describing the coupling of a gluon to a scalar field, whereas the other two terms originate from gluon 
spin or magnetic moment. ^plj,^iq,Pi,P2) coincides with the BFMFG bare vertex involving one background (g) and 
two quantum {pi,p2) gluons |ll|. Evidently the above decomposition assigns a special role to the q-leg, and allows 
^Fa/iiy to satisfy the Ward identity 



<i"^fUM^Pi^P2) = ip'i -Pl)9f.^, (2.4) 

where the right hand-side (RHS) is the difference of two-inverse propagators in the FG, and vanishes "on shell", 
i.e. when pf = P2 = 0. As has been explained in detail in [^jl4|, and as we will discuss extensively later on, the 
splitting of rL°^^(g,_pi,P2) into r^Q^^('7,Pi,P2) and Pp^^^ (9,^1,^2) given in Eq. (2^) has a natural interpretation 
in the context of the tree-level process q{P)q{P') —>■ g{pi) + g{p2) (annihilation channel), leading to an interesting 
connection with the optical theorem. 
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Consider next the S'-matrix element for the quark ((?)-antiquark {q) elastic scattering process q{P)q{P') q{Q)q{Q') 
in QCD; we set q = P' — P = Q' — Q, and s = is the square of the momentum transfer. One could equally well 
study the annihilation channel, in which case s would be the centre-of-mass energy. We will work in the RFG; this 
constitutes no loss of generality, as long as one studies the entire gauge-independent process. It is a straightforward 
but tedious exercise to convince one-self that through pinching, i.e. by simply exploiting the fundamental WI of Eq. 



(2.10), one can arrive at the set of diagrams of the RFG starting from the set of diagrams at any other value of ^ (see 
for example pO| ) , or even from the diagrams corresponding to non-covariant gauge-fixing schemes | ]3l[ | . 
Let us define 

ti^,y{q) = q'^9^il^ - qt^iv , (2.5) 

rf(9)-^, (2.6) 
q^ 

Soip)^-^, (2.7) 
p — m 

Ptiu{v) ^lp^So{p)-ii, . (2.8) 

^A'j^('Z) is the dimensionful transverse tensor, and So{p) is the tree- level quark propagator. In what follows we will 
use the short-hand notation [dk] — with d = 4 — 2e the dimension of space-time and /x the 't Hooft mass. 

Furthermore we define the scalar quantities 

Mq,k)^g^CA[k\k + qf]-\ 
J2{q,k)^g^CA[k\k + qf]-\ 

MqJ, k) = '-g^CA[k^{k + e)\k + e - qf]~' , (2.9) 

where Ca is the Casimir eigenvalue of the adjoint representation ( Ca = N for SU{N)). The quantities J2{q,k) and 

Jsiq, £, k) will be used in later sections. 

We then implement the vertex decomposition of Eq. ( ^.2|) inside the non- Abelian one- loop quark-gluon vertex graph 

of Fig.la, to be denoted by Fq ^'"'"''(Q, Q'), where now pi^ = -/c^, p2i, — [k — q)v The Tpl^^^{q,pi,p2) term triggers 
the elementary Ward identity 

^ = (^+ @-"^)-(0-"^); (2.10) 
thus, a self-energy like piece is generated (Fig.lc), which is to be alloted to the conventional self-energy. In particular, 

Fii)^'-^(g, Q') = r(i)^"'^''(Q, Q') + i v},ZiQh^ 

+x['JiQ, Q'){Q' - m) + (0 - m)X^'J{Q, Q') , (2.11) 



where 



r«''-^(g,Q') = J [dk]J,{q,k)TP^^^{q,-k,k~q)p'^-^m, 
ypLil) = ^ J[dk]Mq,k)g^^, 
X['J{Q,Q') = I [dk]J,{q,k)j^So{Q'), 

Xi'J{Q,Q') = / [dk]J,{q,k)SoiQ'ha ■ (2.12) 



The terms in the second line on the RHS of Eq.(2.11) vanish for on-shell external fermions. The (dimension-less) 
self-energy- like contribution ^ ^p^aa-il)^ together with another such contribution arising from the mirror vertex (not 
shown), after trivial manipulations gives rise to the dimensionful quantity 

^'^pU<l)^V^Zmil)- (2.13) 
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^Pa/sil) will be added to the conventional one- loop two-point function n^^^(g), to give rise to the the PT one-loop 
gluon self-energy ni^^(g): 

nS(9) = nS(g) + nW^(g). (2.14) 

In particular, suppressing color indices throughout, we have that n|^^^(q) is given by the graphs of Fig 2a and Fig 2b, 
namely 

n^l?) [dk]Mq, k)L^p{q, k) , (2.15) 

where 

L^p{q, k) = T^^'^^''{q, k, -k - q)Tf^^{q, k, -k - q) - 2fc„(fc + q)p , (2.16) 
and thus the PT one- loop gluon self-energy u''^^{q) (Fig. 2) assumes the closed form Q 

nSC?) = IJ [dk]Mq, k)L^^{q, k) , (2.17) 

with 

L^p{q, k) = rf^^'iq, k, -k - q)Tfp,p{q, fc, -k - q) ~ 2{2k + q)^{2k + q)p . (2.18) 

Notice that in general both ^^^^{q) and ^''pl^piq) on the RHS of Eg. ( 2. 14 ) depend explicitly on the GFP in such a 
way as to give an GFP-independent sum. 

In addition, gauge-invariance is encoded in the WI 

qaU^^liq) = q^n^^liq) = . (2.19) 

The following important points have been discussed in detail in the literature (i) Hap^q) is independent of the 
gauge-fixing parameter in any gauge-fixing scheme, (ii) As happens in QED for the photon self-energy p2| , the 
gluon self-energy n'^^ (g) captures the leading logarithms of the theory at that order ||l|,^ ; therefore the coefficient 
in front of the single logarithm coming from 11^^^ (g), coincides with the first coefficient of the QCD /? function [ p3[ . 
(iii) ^ap{q) can be Dyson resumed, following the diagrammatic algorithm presented in ||l^. (iv) The combination 

'^e//(9iA*) ^ 5'^(m)^(9/m)7 where A(q/ ii) = [1 — li*-^^ (gZ/i)]^^ is a renormalization-group-invariant quantity, and 
constitutes the non-abelian analogue of the QED concept of an effective charge pP,p7|. Additional properties for 
the gluon two-point function have been presented in the literature for the case of non-abelian gauge theories with 
Higgs mechanism |Q,^|2]-|l^] , and the (non-renormalizable) Kunimasa-Goto-Slavnov-Cornwall ||3^ massive Yang-Mills 
model [3^ ]. 

The PT quark-gluon vertex f'^\Q,Q') is the sum of the non-abelian and abelian one-loop graphs, shown in 
Fig. 3a and Fig. 3b, which we will denote by Pa ^'""''(Q, Q') and Ta^'"^'' {Q, Q'), respectively. In addition to being 
GFP-independent, by virtue of Eq. (O) fi^^((5, Q') satisfies the following QED-like WI 



qc.fi'\Q,Q') = &'\Q) - S(1)(Q'), (2.20) 

where S^^) is the PT one-loop quark self-energy, which coincides with the conventional one computed in the RFG. 
@- 

Finally, the PT one-loop ?i-point functions coincide with those computed in the BFMFG ("tilded" quantities) |l^] 
, i.e. 

nS('?) = nS(?,eQ = i) (2.21) 

fW(Q,Q')=rW(Q,g',CQ = l) (2.22) 
S(i)(Q) = S(i)(Q,eQ = 1) - ^^'HQ,C = 1) (2.23) 

In addition, exactly analogous properties have been established for the one-loop gluon three-point function |^ and 
four-point function ( p6| and second paper in pf). 
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III. THE FULL TWO-LOOP CONSTRUCTION 



Here we present the full two-loop construction. The basic observation is the following: if one carries out the 
decomposition for the bare three-gluon vertex described in Eq. (|]^) to all external vertices (to be defined in sub- 
section IIIB) appearing in the Feynman diagrams contributing to the two-loop S'-matrix element for the process 
qq — > qq, then two-loop sub-amplitudes will emerge, with precisely the same properties as the one-loop PT effective 
Green's functions. 

Throughout this section we have used the following formulas, valid in dimensional regularization: 

0, 



J [dk]{:2k + q)^,h{q,k) 



\-n^ (k'^\ 

[dfc]— ^=0 iV = 0,l,2,... (3.1) 

The last relation guarantees the absence of tadpole and seagull contributions order by order in perturbation theory. 
In the two-loop calculation presented in this paper only the case = 1 is relevant. Notice however that nowhere 
have we used the slightly subtler dimensional regularization result 

which is often employed in the literature. We also use the group theoretical identities 

jaex jbex (Jj^S^^ 

jaxra jbmn jcnx j|_ ^^jabc 

where t° are the gauge group generators in the fundamental representation; in the case of QCD r° — A''/2, where A" 
are the Gell-Mann matrices. 



The identity 



-r(O) p(0)/3m'' _ p(0) p(0)/3Mi' , p(0) y.(a)|3^LV -p(O) p(0)/5A"' p(0) p(0)/3M!' 
^ aiiv^ — ^ Fafiy^ F ^ Pafii^^ F Fafiv'^ P ^ Panu^ P 



_ p(0) -p(0)/3Mi/ -p(O) p(0)/3m<' , p(0) p(0)/3M!^ _ p(0) y.(0)l3^iv 

— ^ Fafiv^ F T~ ^ Pativ^ ^ Q^i/^ P ^ Pafii^^ P W-^> 



may also be found useful at intermediate steps. 
The two-loop integration symbol 



will be suppressed throughout. We define the following quantities 



i^^-r I I 7^7^ (3.5) 



ih = 9^C\[e{i - qfk^k -f if(k + i - qf\-^ , 
ih^g^Cl\i\i-qfk\k^qf\-\ 
ih=g'C\\i\i-qfk\k^lfr\ 
j/4 = g''C\\i^e{i - qfk^{k + £)2]-l , 

ih^g''C\{ek\k + qfY\ (3.6) 



which will be used extensively in what follows. 
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A. The one-particle reducible graphs 



As has been explained in detail in the resummability of the one-loop PT self-energy requires the conversion 
of one-particle reducible (IPR) strings of conventional self-energies H^^^ into strings containing PT self-energies 
The process of converting conventional strings into PT strings gives rise to left-overs, which are effectively one-particle 
irreducible (IPI), and must be alloted to the genuine IPI two-loop structures. Various self-consistency arguments 
supporting the validity of this method have been presented in the literature js^] ; as we will see at the end of the third 
sub-section, the extension of the PT to two loops provides the ultimate test for the self-consistency of this procedure. 
It is straightforward to establish that the set of IPI graphs (Fig. 4a - Fig.4d) may be converted into the equivalent set 
of IPI PT graphs (Fig.4e - Fig.4h), up to some missing pieces: 



(46) + (4c) 



(4a) 



(4e)- 
(4/) 



R 



(2) 
Pa/l 



iAg) + iAh)^Fl?l{Q,Q') 



(3.7) 
(3.8) 



with 



nW(g)T/A^''(g) + ^n«^(9)T/«''(g) 
(g)d(g)r« (g, Q') + r;2j(Q, Q') 



(3.9) 
(3.10) 



with 



r^'j(Q, Q') ^ Q')s(^nQ) + ^L^Q, (Q') ■ 



(3.11) 



The above terms originate from carrying out the vertex decomposition of Eq. ( |2.2D at all conventional IPI diagrams. 
For example, the term Yp^{Q^ Q') (Fig 5b) originates after imposing Eq. ( pT^ ) on the diagram of Fig. 5a ; in addition, 
one obtains the PT counterpart of Fig. 5a, namely Fig.5a^, and the graph of Fig. 5c, which is part of Fig.4f. The terms 
Rp^apiq) has been derived in detail in |lj]. Notice that 



R 



(2) 
Pa/3 



(3.12) 



B. The two-loop vertex 

In this sub-section we will demonstrate the construction of the two- loop PT quark-gluon vertex Fq, (Q, Q'), which 
turns out to have the exact same properties as its one-loop counterpart F^ {Q, Q'). At the same time we will determine 

f2) (2) 

the two-loop propagator-like contributions VpJ^7o-, which will be subsequently converted into Hp^^, i.e. the two- loop 



version of Hp^^ of Eq.( ^.13 ). In addition, out of this procedure the terms Yp^{Q, Q') of Eq.( 3.11 ) will emerge. 

The construction proceeds as follows: The Feynman graphs contributing to Ta {Q,Q') can be classified into two 
sets, (a) those containing an "external" three-gluon vertex i.e. a three-gluon vertex where the momentum q is 
incoming, as shown Fig. 6 and Fig. 7, (b) those which do not have an "external" three-gluon vertex. This latter set 
contains either graphs with no three gluon vertices (abelian-like) , or graphs with three-gluon vertices whose all three 
legs are irrigated by virtual momenta, i.e. q never enters alone into any of the legs; such would be for example the 
abelian graph of Fig. 3b, if one was to insert a one-loop self-energy correction to the internal gluon line [Q). Of 
course, all three-gluon vertices appearing in the computation of the one-loop S'-matrix are external, and so are those 
appearin g in the IPR part of the two- loop S'-matrix (see previous section). Then one carries out the decomposition 
of Eq. ( ^.2|) to the external three-gluon vertex of all graphs belonging to set (a), leaving all their other vertices 
unchanged, and identifies the propagator-like pieces generated at the end of this procedure. 

The calculation is straightforward, but lengthy; it is more economical to identify the sub-structure of the one-loop 
three-gluon vertex T^al.v{q-,Pi,P2) (Fig K) nested inside the two-loop graph, Fig(Nal), and use the results presented 



in section VI, Eq. (6.35). To that end we must set pi — £ and p2 — > £ — g, and J3{q, —pi, k) Jaiq, k), and 
rewrite the Fp'^^^^ of Eq.( 6.35 ) as follows 
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^3 [r(°U^ -q,k,-k~£ + q) + Kg^p]iP - V^Zil) 
h [r|Sp(-^, k + ^, -fc) + k^g^p]{l ~qY-l Vj^'^iq) 



{i - q)u 



h [T^^U^ - q, k, -k-£ + q) + + 1 d{q)q^V^'Uq) 



Thus we have (we omit the external spinors, and use Q' ^ Q + q) 



(6ai) = (6af ) + ln'-^}^\q)d{q)T^'^'''^''{Q, Q') 



2/25a. - h [r^l i-k,-£,k + £) + k^g^p] {£ - qY 



where 



(65i) 
(66i)' 
(6ci) 
{^n) 



D3 



+Di +D2 + D3 + Di 

(66f ) + ^hgaal" -Di+D5 + De 

(6cf) + i/4L„a(^,fc)7" 
(7rf) + in«^(g)d(g)r(^)-''(Q,Q') 



Da 



- fi 

D,, + De + Dj + Ds 



-hK>M, k, -£-k) + k,g^^]{YS{Q + £ + fc)7'^} , 
h[r^°U^, ^£ -k,k) + k^g^,]{j-'S{Q' - £ - fc)7^} , 

-lhqa{l^S{Q + £)r}, 
1 



D4 = 7;hqa{l^.S{Q' + ihn 

D5 



2 -oiiJL L ' M*^ 

-hilaSiQ - kY^S{Q -k- £)r} , 
De = h{laS{Q - k)^^S{Q + £)r} , 
D7 = -h{lf.S{Q -k- £)rS{Q - fc)7c} , 
Ds = h{l^S{Q + £h^S{Q - fc)7a} , 



(3.13) 



(3.14) 
(3.15) 

(3.16) 

(3.17) 

(3.18) 



(3.19) 



Before we proceed the following comments are warranted; 

(i) In the above formulas appropriate shiftings and relabellings of the integration momenta have been carried out, 
in order for the answer to be expressed in terms of the five basic denominators li defined in Eq.( |3.6| ). 

(ii) The topologies of Di, D2, D3, D4, D5, and Dg are shown in Fig. 604, Fig. 605, Fig. Gag, Fig. 607, Fig. 664, and 
Fig. 665, respectively. (65i)' corresponds to the figure obtained from Fig. 661 by drawing the internal three-gluon 
vertex on the other leg of the external three-gluon vertex, i.e. the leg which hooks onto the external spinor carrying 
momentum Q' (not shown); D7 (shown in Fig.T/i) and £)§ (not shown) arc the analogues of D5 and Dq for the (661)' 
topology. Furthermore, the second and third term on the RHS of Eq.(3.14) are depicted in Fig. 602 and Fig.Gas 
respectively, the second term on the RHS of Eq.(3.15) is shown in Fig. 662, the second term on the RHS of Eq.(3.17) 
is shown in Fig.6c2, and the second term on the RHS of Eq.(3.18) in Fig.7ca. 

(iii) The graph in Fig.7r2 is accompanied by the graph with the abelian vertex correction on the other side (not 
shown). The graph containing the bare four-gluon vertex is not shown either. 
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(iv) Notice that the propagator-hke part in Fig. 7a is connected to the rest of the graph by a factor gap] since this 
term will become part of the IPR graph of Fig. 4c it must be supplemented a longitudinal component qaqp- This 
contribution is proportional to the terms shown in Fig. Tag and Fig. Toy (remember that the latter are proportional to 
get), but notice that the color factors are different; the terms shown in Fig. Gag and Fig. 607 are proportional to C\, 
whereas the term missing in Fig.(Sa) will have both C\ and CaC/, where C/ is the Casimir for the representation 
of the (external) quarks. It is easy to show that these terms vanish when the external fermions are on shell, which is 
our case. In order to prove that one does not have to assume that qa hits a conserved current on the other side of the 
graph; instead one notices that the terms in question are proportional to [I]^^-'(j^,m) — T,^^\Q' ,m)], which vanishes 
when Q —Q' = ni. 

(v) Notice the appearance of propagator-like terms (Fig.6a3, Fig. 662, and Fig.6c2) 
Adding Eq.( |3.14 ) - Eq.( |3.18D , by parts, and using Eq.( f3.11[ ), we find 



(2) 



(Q,Q') + 2^pi7^ 



with 



VplM = k) + (2/2 + h)gc.a - h 



kagap + TpCTQ ( — fc, — ^, k + £) 



(3.20) 



(3.21) 



The interpretation of the three terms appearing on the RHS of Eq.(3.2C) is as follows: 

(i) The term \F^]^{Q,Q') is half of the vertex-like part necessary to cancel the corresponding term appearing in 
Eq.(3.8), during the conversion of conventional IPR graphs into their PT counterparts. The other half will come from 



the mirror vertex (not shown). 

(ii) ^V^QCT^"^ ^^"^ total propagator-like term originating from the two-loop quark-gluon vertex; together with the 
equal contribution from the mirror set of two-loop vertex graphs (not shown) will give rise to the self-energy term 



n 



(2) 

Paf3 



(3.22) 



which will be part of the effective two-loop PT gluon self-energy, to be constructed in the next sub-section 

(iii) Ta {Q,Q') is the PT two- loop quark-gluon vertex; it coincides with the corresponding two- loop quark-gluon 
vertex computed in the BFMFG, i.e. 



ri'HQ,Q') = r(2)(g,Q',eQ = i) 



(3.23) 



as happens in the one-loop case , Eq.( ^.22 ). Either by virtue of the above equality and the formal properties of the 
BFM, or by means of an explicit diagrammatic calculation where one acts with qa on individual diagrams, one can 
establish that Fq (Q, Q') satisfies the following QED-hkc WI 



qc.f'i\Q,Q')^&'\Q)-t^'HQ') 



(3.24) 



which is the exact two-loop analogue of Eq.( 2.20| ). I](^^((5) is the two-loop PT fermion self-energy which satisfies 



£(2) (g) ^ (g, = 1) = s^'^ (g, e = 1) 



(3.25) 



Again, this is the precise generalization of the one- loop result of Eq.( 2.23| ). At this point this result comes as no 
surprise, since all three gluon vertices appearing in the Feynman graphs contributing to E^^-'(g,^) are internal; 
therefore, at ^ = 1 there will be no pinching 

We emphasize that the above result is non-trivial; indeed, even if one accepts that the appearance of the first and 
third term in Eq.(3.20), for example, could be forced, there is no a-priory reason why the remainder should turn out 
to be purely self-energy-like. Notice also that in deriving the above results no integrations (or sub-integrations) over 
virtual momenta have been carried out. 



C. The two-loop self-energy 



^ (2) (2) 
The construction of n]^^(q) proceeds as follows: To the conventional two-loop gluon self-energy n^^(q) we add two 



(2) 

additional terms; (i) the propagator-like term np^^(q) derived in the previous sub-section, Eq. (|3.22| ), and (ii) the 
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(2) I I 

propagator- like part ~Rp anil) given in Eq.(3.£), stemming from the conversion of the conventional string into a PT 



string; this term must be removed from the IPR reducible set and be alloted to Tl^^^^{q)^ as described in ||T2|. Thus, 
na/3(9) reads 



(3.26) 



Up to minor notational modifications, this last equation is in fact identical to Eq.(3.5) in the second paper of p2[ 
except that now we know the exact closed expression for the term ^paaio)- 



It is a lengthy but relatively straightforward exercise to establish that in fact 

nS('/) = nS(9,eQ = i) 



(3.27) 



To see that in detail, we simply start out with the diagrams contributing to nj^^(g) and convert them into the 

— (2) 

corresponding diagrams contributing to Il^^{q, = 1); in doing so we only need to carry out algebraic manipulations 
in the numerators of individual Feynman diagrams, and the judicious use of the identity of Eq.(3.4). The individual 
diagrams yield: 



(8a) = (9a) + (9m) + (97i) + (9o) + {9p) + (9g) + (9r) + (9s) 



V'fc^rW (g, k + e-q,-k-e) + ii-q-£^)ik + i- q)^{-£ - q)p 



+ r(l(-fc, -I, k + em - q)%iq) - (2/2 + i3)tM 
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2 ^ 



lo.{l-q)0 + {i-q~f)go.p 



+l2[Lc.fi{q, k) + 2ka{k + q)p] - hiLcnii, k) + 2fc„(fc + C)fi\ 
21 

3^a^/3 - -^qa^l3 + qaqi3 



{8b) + (8c) = (96) + (9c) + -h 

(8d)-(9d) + ^/2ia/3(g) 
o 



e4i-q)0 + ie-q-e^)gc.^ 



3) + (8/) = (9e) + (9/) - /i V'fc'rL") (g, k + e-q,-k-i) 



2ka{k + q)fj + -qaqf3 



2ka{k+i)i3 + -qa 



(85) + i8h) - (9.g) + (9/1) + hL^pie, k) - /4i„.(^, k)t%) 

(8z) = (9i) -h{l-q- e){k + I - q)^{^-e - q)p ~ ^hq^qp + \hqaip 

(8j) + (8fc) = (9j) + (9A:) + /sp^a^/j - S^ag^j + g^g^] 
(8^) = {U) 



(3.28) 



Adding the above equations by parts we find 



ni'i(a) = n^1(g,eQ = i) + /2 



;(2) 



Lap{q,k) + 3tap{q) 



hL^,{i,k)tll{q) 



fc,5ap + r(l(-fc,-^,fc + €) 



{l-qYt{q)^(2h+hMq) 



(3.29) 



Using Eq.(3.21) and Eq.(3.12), Eq.(3.29) becomes 
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nS(9) = ni^j(9, eg = 1) - n(?L,(9) + RpUi)- (3-30) 



From Eq.(^30|) and Eq.(p^) we arrive immediately to Eq.( ^.27| ). Again, no integrations over virtual momenta need 



be carried out, except for identifying vanishing contributions by means of the formulas listed in Eq.(3.1). 

Since we have used dimensional regularization throughout and no integrations have been performed the results 
of this section do not depend on the value d of the space-time; in particular they are valid for d — 3, which is of 
additional field-theoretical interest pO[ . Clearly, when c? — > 4 the renormalization program needs be carried out; this 
will be the subject of the next section. 

IV. RENORMALIZATION 

In this section we will carry out the renormalization for the two-loop PT Green's functions constructed in the 
previous section. There is of course no doubt that if one supplies the correct counterterms within the conventional 
formulation, the total S'-matrix will continue being renormalized, even after the PT rearrangement of the (unrenormal- 
ized) two-loop Feynman graphs. The point of this section is to show a stronger version of renormalizability, i.e. that 
the new Green's function constructed through the PT rearrangement are individually renormalizable. The general 
methodology is as follows: We start out with the counterterms which are necessary to renormalize individually the 
conventional Green's functions contributing to the two-loop S'-matrix. Then we will show that by simply rearranging 
them, following the PT rules, we will arrive at renormalized two-loop PT Green's functions. 

We will use the following notation: Zi is the vertex renormalization constant for the conventional quark-gluon vertex 
Ta, Zi is the wave- function renormalization for the (external) quarks, Za the gluon wave- function renormalization 
corresponding to the conventional gluon self-energy 11, Z a the gluon wave-function renormalization corresponding to 
the PT gluon self-energy 11, Zj, is the vertex renormalization constant for the conventional one-loop three-gluon vertex 
Ta^i^, and Z-^p the vertex renormalization constant for the P^lj^^, Z^ is the ghost wave- function renormalization, and 
Z\ the ghost-gluon vertex renormalization constant. Equivalently, one can carry out the renormalization program 
using appropriately defined counter-terms. The corresponding counterterms, which, when added to the above n-loop 
quantities render them UV finite, are, respectively k'^^ , K^^^ , k'^^^ , -f^^"'' , ^3"' , ^3^-^ -^2""* ' ^"^^ . In addition, 
mass counterterms 8m must be supplied if the quarks are considered to be massive. In a moment we will also introduce 
the counterterm i4rjr\ which renders ultra-violet finite. Notice that because of the QED-like WI it satisfies 
becomes ultraviolet finite when the counterterm K2 is added to it. The Z& and the Xs are related as follows : 

= l + ^Xp') z = 1,2,3 



Za = \^Y.^_ 
Z^ = 1 + ^X 



A 



U) 

A 



= l + ^Xp') z = l,2 

(5m = ^(5m(^') (4.1) 
i=i 

We first begin with the IPR part of the S'-matrix. It is more convenient to work with dimension-less quantities; to 
that end we define the dimension-less gluon self-energy II'^' simply through = Tl^^hap. In order to renormalize 
the extra pieces which must be alloted to the conventional IPR graphs in order to convert them into their PT form. 



i.e. the iRp + Fp terms in Eq. (B.9) - Eq. (3.11), we must have (we use the supescript "R" to denote renormalizaed 
quantities) 



4' 



11 



iR 



(2) 



F. 



(2) 



(2) 



u. 



(2) 



(4.2) 



with 



(2) 



U. 



V, 



(1) 



K. 



(1) 



if 



(1) 



(2) 



K 



(1) 



(1) 



K 



(1) 



n 



(1) 



(1) 



-(1) 



V, 



(1) 



,(1) 



(4.3) 



Thus, the terms U^^"^ and U2^' need be supphed. Notice that the terms contained in t/j"^'' are momentum independent 

(2) 

(quadratic in the counterterms) , whereas those contained in Ul or momentum-dependent (hnear in the counterterms) . 
The latter will cancel against parts of the one-loop counterterms appearing inside the two-loop expressions for the 
conventional self-energy (Fig. 10) and vertex (Fig. 12), cancelling their sub-divergences. As for the former, they will 
become part of the final renormalization counterterm of the two-loop PT self-energy, i.e. the counterterm necessary 
for cancelling the remaining divergence, after the sub-divergences have been taken care of. Another way of saying 
this is that, since the extra terms Rp and Fp will be alloted to the PT gluon-self-energy and vertex, so should the 
counter-terms necessary to renormalize them. 

(2) 

We next show how the terms in [/, are to be accounted for, and, at the same time, derive some useful relations 



r(2) 



(2) 



among the various counterterms. The QED-like WI of equations Eq. ( 2.20| ) and Eq. ( ^.24 ) relating Pq^ and and 
Ta^ and S^^-*, respectively, imposes the following QED-like relation between the renormalization constants Zi and 
Z2, up to order g^: 



Zi — Z2 ■ 



In addition, from the WI of Eq. ( 3.36 ), we have to order (at least) 

Z3F = Za ■ 



(4.4) 



(4.5) 



Notice also that Eq. ( 6. 36] ) dictates that the ultraviolet-divergent part of Pi?^^^ is proportional to ^p^^^^, rather than 



I^Fa^ii/' ^1 been the other way around there would be no longitudinal ultraviolet-divergent pieces on the RHS 
of Eq. ( 6.36 ). As we will see, this "mismatch" will generate the pieces which in the BFM language give rise to the 
gauge-fixing renormalization of the vertices (Fig. 10a and Fig. 10b) 

Furthermore, the renormalization constants before and after the PT rearrangements are related to the gauge 
coupling renormalization as follows: 

'72 r72 ry — 2 ry — 1 

— ^1^2 '^A 



Z\Z2 '^Zj^ 

z-.\ 



(4.6) 



where Eq. (4.4) has been used. After substituting the expressions given in Eq. ( [4.l| ) into Eq. ( |4.6| ) and equating 
powers of the coupling constant g we arrive at the following relations for the corresponding counterterms 



'(1) 



(1) 



-K. 



K 



(2) _ 7^(2) 



'(2)x 



'(1) 



K. 



(1)1 



Substituting the relations 



into the above equations we obtain 



K 



K. 



1,2 



K 



(1) _ 7^(1) 



^A 



(1) 



K 



A 
(1) 



K 



(2) 



u. 



(2) 



(4.7) 
(4.8) 

(4.9) 

(4.10) 
(4.11) 
(4.12) 
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Notice that if Z\ = Z2, or equivalently, Kp^ = 0, then Za = Za, which is simply the QED case. 

(2) 

It is now clear what the role of the U2 terms is; they must be added to the conventional two-loop self-energy 

(2) . (2) 

coimterterm Kj^ and the part of the vertex counterterm corresponding to Vp ; these two counterterms are already 

(2) 

present, whereas the term U2 must be borrowed from some other part of the S'-matrix, in order for the initial 

equation Eq. (O) to be enforced. 

(2) 

We next turn to the terms. We will show that they will cancel precisely against terms originating from the 
rearrangement of the graphs shown in Fig. 10 and Fig. 12, in order to convert them into the graphs shown in Fig. 11 
and Fig. 13. The former set contains the counterterms necessary for cancelling the one- loop sub-divergences inside the 
conventional two-loop gluon self-energy and quark-gluon vertex, the latter the counterterms needed for the two-loop 
FT gluon self-energy and quark-gluon vertex. 

We begin with two-loop gluon self-energy shown in Fig. 10. From the Slavnov- Taylor identity we have that 



Z3 _ Zi 
Za Z2 



from which we obtain |43l 



This relation is important for what follows. 
Then we have for the graphs of Fig. 10: 



(10a) + (106) ^ Ki'^ I [dk].h [q, fc)r(%r('2, 

(10c) = (lOci) + (IOC2) 



with 



(4.13) 



= (4.14) 



(lOd) + (lOe) = k[^^ J [dk]Ji {q, k)k„{k + q)p 

(10/) + (IO5) = -X^'^ [[dk].h{q,k)k^{k + q)p (4.15) 



(lOci) ^ -/i « j [dk]j,{q, fc)rW,rJ2, 

(IOC2) = K^}^ j [dk]J2{q, k)k^KT^°lTf^^ ■ (4.16) 
J2 has been defined in Eq. (^.g|). Thus, using Eq. (4.14) and Eq. (2.15) we obtain 

[(lOa) + (106) + (lOci)] + [(lOd) -f (10/)] + [(lOe) -f (lOg)] = K^p^H^^liq) (4.17) 
Next, we write (IOC2) as follows 



(IOC2) = (11a) -f (116) -t- (He) - 2i„^(9) j [dk]J2{q, k)kf3k'' (4.18) 

where 

(11a) + (116) = xi^) I [dk]Mq, k)rf^,,r%'"' 

(11c) = Xi^) I [dk]J2{q, k)k^k,Tf^^yZ'''' (4.19) 

(2) 

We next convert the renormalization constants for the conventional two-loop quark-gluon vertex Fq (Fig. 12) into 
those necessary for the PT 2-loop quark-gluon vertex Fq finite (Fig. 13). We have: 
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'(1) 



^(1)t/(1) 



(13c) + (13d) 



ii^«)r« = (i3a) + ix«r« 



(12a) 

(12c) + (12d) = K\''V, 
(12e) 

(126) - -K^^^ I [dk]Mq,k)Tl^ltl{k)p^''iQ + k) 



1 



(1)t/(1) 



(13e) 



(136) - K^]^^ J [dk]J2{q,k)t^^{k)-f^ 

(136) - ^i^i'V^'^ + K^to^Aq) f[dk]J2{q, k)kpk''^l^ 



(0) 



(4.20) 



Notice that the last term on the RHS of Eq. ( 4.20| ) , together with an equal term from the mirror vertex graphs (not 
shown) will cancel against the last term in Eq.( 4.18| ). Using that k[^^ — ^-^2^' — ^i-^p'' + -^2^'') finally arrive at 



^(Fig.lO) + 2if(_Fig.i2) — ^(Fjg.ll) + 2iir(Fig.l3) + ^^1 



(1) 



(4.21) 



where the factors of 2 multiplying K^pig 12) and K^pig i^) account for the mirror vertex contributions. Evidently, the 

PT rearrangement gives rise to a term as announced. 

The counterterms resulting from the above rearrangement are exactly those required to cancel all sub-divergences 



inside fa''(Q,(5'); the latter coincide of course with the sub-divergences inside Ta' {Q,Q' ~ 1). To demonstrate 
that the counterterms shown in Fig. 13 are in fact identical to those obtained when carrying out the BFM renormal- 
ization program as explained in i.e. renormalizing only the background gluons, the external quarks, the coupling 
constant g, and the GFP (^q), one may proceed as follows: (i) start with graph (13b), and separate the tcr^{k) into 
the part proportional to and the part proportional to knk^\ the second part is simply the gauge- fixing renor- 
malization to the self-energy, as explained in (ii) Half of the piece proportional to must be added to (13a); 



(2), 



the latter is proportional to r^"^ Using that ^^3^ = K)^' (from Eq. (4.5)), the total contribution is proportional 



to Kj^ {T'^'^^ — Tp ) ~ Kj^ Tp \ which is the contribution from the gauge-fixing renormalization of the elementary 
BFM three-gluon vertex, (iii) The remaining half proportional to g^p, from step (i) must be split equally between 
(13c) -I- ^(13e) and (13(i) -f i(13e). Each of these two combinations will then give a contribution proportional to 



K^i^ - 1^2' - hK'x' = -\K^A + 1^2' = Kl'' + \K)^'', where Eq. (|4J|), Eq. ( [4.10D , and Eq. (f4.12| ) have been 
used in the last steps. The wave-function renormalization for the external fermions will then cancel ^K^\ and k\^^ 
will be re-absorbed into the gauge-coupling renormalization. 

Notice that the correspondence between the PT and the BFMFG Green's functions established in the previous 
section persists after renormalization; the resulting expressions are the BFMFG renormalized quantities, as derived 
in ||]. An immediate consequence of this is that the coefficient multiplying the logarithm of the PT two- loop 
self-energy is equal the the two-loop coefficient of the QCD [3 function i.e. (34/3) C\. It is also interesting to 

see how the PT rearrangement leads into the interpretation of counterterms generated from the QCD Lagrangian in 
the renormalized Feynman-gauge |42|, as GFP renormalizations of the BFM Lagrangian B . 



'(1) 



V. THE ABSORPTIVE CONSTRUCTION: GENERAL FORMALISM 

In the next two sections we will show in detail how one may construct the two-loop PT effective Green's functions 
using unitarity and analyticity arguments. This derivation generalizes the method first presented in fl^ and |lj] 
for the one-loop case, and constitutes a non-trivial self-consistency check for the entire approach. In this section we 
will set up the formalism, and discuss in detail the one-loop case, which will serve as the general paradigm for the 
two-loop generalization, to be presented in the following section. Apart from some minor modifications, in the next 
two sections we will adopt the notation used in |p^ . 

The optical theorem for the case of forward scattering assumes the form 
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^m{a\T\a) = 1 ^(27r)V4)(p, -p,)(*|T|a)*(»|T|a) , (5.1) 

i 

where the sum should be understood to be over the entire phase space of all allowed on-shell intermediate states 
i. After expanding the T matrix in powers of g, i.e. T = X]n=2 T^^\ we have that 



In particular, if in the initial states we have a q-q pair, i.e. \a) = \qq) we have for the first non-trivial orders, n = 4 
and n = 6 

'^m{qq\Ti%q) = ^( |) / {dPSU{2g\T^%q)* {2g\T^%q) , (5.3) 



and 



^m{qq\T^%q) = ^ (^^ j / idPSh,{3g\T^%q)* {3g\T^%q) 

+ ^(^) / idPS)2,2^e(^{2g\T^%q)*{2g\T^%q)^ , (5.4) 

respectively Q]. {dPS)2g and {dPS)3g stand for the two- and three-body phase space for mass-less gluons, respec- 
tively. Next we introduce the short-hand notation 

=3m(gg|T["l|gg), 

Tjii:^ = {mg\T^%q) m,fc = 2,3,... (5.5) 



Then, Eq. ( p.3[ ) and Eq. (5^) become respectively 



(5.6) 



1 / 1 



2 \ 3! 



-4" +4''- (5.7) 
The quantities defined above have the explicit form 

= I (^) K'f + 4'tt^'^'^'(pi)^"^'''(p2)[rif + r|i];:t , (5.8) 

4" = I (^j^ K'f + T^X%P'''''ipi)P''-''iP2)P'''''iP3)[r^f + , (5.9) 

4' = (^)^^(k^' + riftt^''^'(pi)^''^'(p2)[rw +ri^^^ , (5.10) 

where we have suppressed the phase space integrations. P^^ is the polarization tensor for mass-less gluons, 

Ptiiy{P,n,r]) = -5/^1. H r]- — (5.11) 

np (np) 

with is an arbitrary four-vector, and 77 a gauge parameter. 
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We will now study the process q{ki)q{k2) — > g{pi)g{p2) at tree-level, using the equations derived above. This study 
will further elucidate the role of the Eq. (2.2) in enforcing perturbative unitarity at the level of individual Green's 
functions, and will set up the stage for the two-loop generalization. 

The BRS symmetry of the original Lagrangian leads to the following identities |Q : 

P2p1T,1'. = 0- (5.12) 

If we split the amplitude into an s-channel and t-channel contribution (s ~ — (ki -f fc2)^ — ipi +^2)^, and 
t — (fci — pi)^ — (^2 — p2)^), the first of the identities in Eq. (5.12) becomes 

(t^s + T^t^ - (^S^,l'^ + " . (5.13) 

The remaining two identities are exactly analogous and will be suppressed throughout. 

The PT rearrangement of the amplitude amounts to a special choice for 72 ^ and 72t , which will be denoted by 
and T^f", respectively. After defining these "Feynman" amplitudes, Eq. ( 5.13| ) reads 

p'.{ri + r-y'^=^sl^^'Us^,^''^yP2., (5.14) 



-lk]F ^ ^[k]F \"'^ _ f Qlk]F{l2} . c.[fc]_F{12}^ '''' 



and to order fc, 

p't(T^r+r^r) =i^s^,T'''' + s^,T'''') P2.- (5.15) 

We will next study the case fc 2 (tree-level, Fig. 14). 

For "on-shell" gluons, i.e. p^ — 0, Pfj.,^{p,n,ri) of Eq. (5.11) satisfies the transversality condition p^P^v — 0. 
Thus one may immediately eliminate the rpQ^jy('Z,Pi,P2) part of Tal,v{q,Pi,P2), which vanishes when contracted 

with the term Pfi^xi{pi)Pvv'{p2), and effectively replace ri°jj.(g,pi,p2) by r^_^^j^(g,pi,p2), as in Fig.l4a. One then 
proceeds by recognizing that the longitudinal parts of the P^^'{pi) and Pvu'(p2) trigger a fundamental cancellation 
p3| , p^ involving the s- and t- channel graphs (Fig. 14a and Fig.l4d), which is a consequence of the underlying BRS 
symmetry. In particular, the action of or P2 on F^^^^^ gives 

Pl^'Pa^uil^Pl^P-i) = iP2 -Pl)aP2y + (Pl ~ pl)9ay + taM) 

P2^Flf,M^Pl^P2) = iP2 -Pl)aPlt, + {pl ~pl)gaf, ~ taf,{q) (5.16) 

The first term on the RHS of cither equation cancels against an analogous contribution from the i-channel graph, 
whereas the second terms vanish for on-shell gluons. Finally, the terms proportional to p2i, and pi^ (Fig. 14b) are 
such that all dependence on the unphysical four-vector and the parameter rj vanishes, as it should. In addition, 
a residual (s-dependent) contribution emerges from these latter terms, which must be added to the parts stemming 
from the g^^ig^y' part of the calculation. In particular ]l^, we have that 

- 5{7™}/™'^'rf(g)ri"i:(9,Pi,P2) , 

r\^]F _ T-[2] 
■^2t ~ ■'2t ' 
C,[2]F{12} _ <-,[2]{12} _ ^ 

^'ll-'2s - W2s ) P2i^ + [J^2 )v , 

rf(4'"')^t = -(A^'l^)^^ (5.17) 

with 
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(A, 



[2]F^ab 
u 
ab 



.[2]F{12} ^ _ 



gmr'diqKiq) 



(5.18) 



and {7°} = igv{k2)T°-^au(ki). 
From the above results follows 



ab / \ ab 

12]\ ^ I gl2]F{12} 



P2v , 



(5.19) 



and thus 



,[4] 



-'2s + ■'2* 



T-[2]F , T-[2 



2(^[2]F{12} 



[2]F{12} 



(5.20) 



with 



" 2 U! 



A 



2 ^2! 
1 / 1 



2 V 2! 



^[2]_F^[2]_F* 
-'2s -'2s 



2 S. 



;2]_F{12} 



^[2]_F^[2]* ^[2]^[2]_F* 
-'2s -'2* I -'2t -'2s 



2]t-[2]' 



s. 



:2]_F{12} 



(5.21) 



In the last step we have defined self-energy (S), vertex (V) and box (B) -like amplitudes, according to their dependence 
on the Mandelstam variables s and t as in the case of a scalar filed theory, or QED; A^g^ depends only on s, Ay^ on s 
and t, and only on t. 

The next step is to identify these sub-amplitudes as the imaginary parts of the effective one-loop self-energy, vertex, 
and box, under construction. For example, for the effective self-energy 11"'' (5) we will proceed as follows: first write 
A^g^ in the form 



then identify 



From the last equation follows 



9mn«(Q)=4i^(g). 



(5.22) 



(5.23) 



^Fl,,uil,Pl,P2)Tplf''{q,pi,P2) - 2(p2 -Pl)a{p2 - Pi) P 



(5.24) 



This last equation leads to a well-defined definition of II^^), without having to resort to an intermediate one-loop 
diagrammatic interpretation: after the two-body phase space integrations has been carried out using standard results 
| p3[ , the real parts may be reconstructed by means of a once-subtracted dispersion relation. Thus, in the absence of a 
full (dispersive) one-loop construction, the li*-^^ so generated does not necessarily have to correspond to the imaginary 
parts of a precisely identifiable set of one-loop Feynman diagrams. Of course, given that the full one-loop construction 
has been carried out, we know that this is actually the case. Thus, when the one-loop (or n-loop) construction for 11 
exists, Eq.(5.22) reduces into a non-trivial self-consistency check. In particular, one must verify that 



C2{Tli%,^Q^l) 



(5.25) 
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where C„{...} is the operator which carries out the n-particle Cutkosky cuts (for n = 2 we have two-gluon and 
two-ghost cuts) to the quantity appearing inside the curly bracket (Fig. 15). This is indeed the case |0, as can be 



directly verified from Eq.(2.17). In fact, the residual contributions originating from the terms proportional to and 
P2 mentioned above correspond precisely to the Cutkosky cuts of the one-loop ghost diagrams (Fig. 15b). 
To fully appreciate the subtlety of the above construction the following comments are now in order: 
(i) Of course, the BRS-driven cancellation for the process qq — > gg takes place regardless of the FT rearrangement 
of the amplitude in general, and the FT decomposition of T^ali' in particular. Indeed, if we had not eliminated 
Fp^,,^(g,pi,P2) but had kept instead the full vertex Fi')l^((?,pi,p2) (as is usually done), the WI analogous to Eq. 



(5.16) would be simply 

P2^^°lM,Pi,P2) = WiPi) ~ t^M ■ (5.26) 



The parts which participate in the BRS cancellation, namely the parts proportional to tau{q) and ta^{q) are thus 
unaffected. What changes after the FT rearrangement is the resulting absorptive part of the effective n-point func- 
tions under construction. So, if one was to define the absorptive part of an effective self-energy keeping the full 
^a^iyiq,Pi,P2), but still exploiting the BRS cancellation in order to elliminate the longitudial terms, one would arrive 
at the imaginary part of the conventional self-energy in the REG, '<smllai3{q); the latter, for one thing, does not 
capture the running of the QCD coupling. 

(ii) One could define the absorptive part of an effective self-energy before carrying out the BRS cancellation. In 
that case one would be led to the absorptive parts of the gluon self-energy in the light-cone gauges; in particular, 
the final answer would depend explicitly on the unphysical quantities n and 77. A dispersion relation would give 
rise to a pathological quantity, since the gluon self-energy computed within the axial gauges is not multiplicatively 
renormalizable, due to its dependence on Ua | |47| . Furthermore, spurious infrared divergences appear in the Feynman 
parameter integrations, which are artifacts and cancel out only when full physical quantities are computed |^^ . 

Thus, one has to first carry out the FT rearrangement at the level of the ^-matrix element, then enforce the BRS 
cancellation at the level of the cross-section, and, only after these two steps, one should define self-energy /vertex/box 
absorptive parts, as one would for a scalar field theory. 

Having set up the formalism and discussed the general methodology, we next proceed with the two-loop absorptive 
construction. 



VI. THE TWO-LOOP ABSORPTIVE CONSTRUCTION 



In this section we will show how to extend the methodology established in the previous section to the two-loop 
case. This construction involves two parts: the first part is the study of the one-loop amplitude for the process 
q{ki)q{k2) 9{pi)g{P2)'i the second is the study of the tree-level process q{ki)q{k2) 9{pi)g{P2)g{pz) ■ As we will 
see in detail the FT rearrangement (at the level of the S'- matrix), will give rise (at the level of the cross-section) to 
the correct Cutkosky cuts. 



A. The one-loop version of F^' and Fp' 

In this sub-section we will first study the conventional one-loop three-gluon vertex, and will see how one can 
arrive at the one- loop generalization of the tree- level vertices F^^^^^ and Fp^^,^, defined in Eq.(^.3D, to be denoted 

by Fp^^j^ and Tp^^^^, respectively. Then we will see in detail why casting the one- loop 5- matrix for the process 
q{ki)q{k2) g{pi)g{p2) into its FT form is crucial for enforcing the optical theorem at the level of individual two- 
loop FT Green's functions. 

We start with t he one-loop 5-matrix element for q{ki)q{k2) —>■ g{pi)g{p2)- 



Then Eq. (|5.15D yields (fc = 4) 



/if 



+ = 4t'™ + 4t'™) P.., (6.27) 
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where 7^^^'^ and T^f^^ are the "Feynman" versions of 7^'^' and 72t; their exact form will be specified shortly. The last 



two equations in ( 5.17 ) become 



Pl\J-2t )nu-{^2t ) P2u - [J^2 )u ■ 



Using Eg. ( 5. 19 ), we have that 



(6.28) 



^2s ' -'2* 



n-mF ^[2] 



-2 5. 



.[4]F{12} ^ ^[4]F{12} 



S. 



[2]F{21} 



(6.29) 



with 



•^25 



^2V 



A. 



[6] 
2B 



2! 



3?e 



-[4]F^[2]F* 
^2s -'2s 



^[4]F^[2] 



3?e 



T-[4]F^[2]' 
■^2t -'24 



J4]F{12} 
'2s 



^[2]F^[4]F* 
-'2s -'2* 



,[2]F{21} 
'2s 



2 5. 



;2]F{12} 



s. 



[4]F{21} 



(6.30) 



From these last quantities one could define, for example, the quantity A!2Sapil) exactly as in Eq. (5.22), i.e 



A^2l-^{r}diq)A':Lpil)diq)m 



[6] 



(6.31) 



The A^2Sapil) must then be such that 



K2) 



(2) 



(6.32) 



The question is what is the correct form of 7^^^' and 7^^, and the corresponding 8^2}^'^'^'^^ and 8^2}^^'^'^^ ; the latter 
quantities are automatically determined once the former have been specified. In particular, if such a rearrangement 
exist, does it correspond to a structure already known from the one-loop PT analysis ? The natural candidate for 
this is clearly the PT rearranged one-loop matrix element for q{ki)q{k2) —>■ g{pi)g{p2), shown in Fig. 16 (individual 
Feynman graphs are shown in Fig. 17, Fig. 18, and Fig. 19); if that were the case, one would begin to discern an iterative 
pattern. As we will see in detail, this is indeed what happens. 

Throughout this section we will suppress the one- loop integration symbol J [dk\ , and will use (see Eq. ( p.9[ ) ) 



J3 = Uq, -pi,k) - -g^CA[k\k - Pif{k + p2Y]-^ 



(6.33) 



We focus on the part of the process involving the conventional one-loop three-glu on v ertex Fajiy ((7,^1,^2), which is 
diagrammatically shown in Fig. 17 . Let us now carry out the decomposition of Eq. (2_^) to the elementary three-gluon 



vertex rS]p((7, k + p2^ —k+pi) appearing in the diagrams of Fig. 17a and Fig.l7d; the parts stemming from the Fp^ 
will propagate towards the remaining elementary three-and four gluon vertices, and will trigger further WI of the 



(0) 



type shown in Eq. ( 5.26 ). Reorganizing the terms thusly generated, one can show that the T^al,v{qTPi,P2) of Fig. 17 
can be written in the form 



^aluiq^Pi^P^) 



d{q)^'PI!M^pi,P2) + r^^Lp^('z,pi,P2) + r^'a^,.(g,pi,P2) , 



.(1) 



(6.34) 



with 
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(6.35) 



The first term is precisely half of the pinch contribution needed for converting Yl^aliq) into Yi^ai/{q), as shown in 
Fig. 16a; the other half will come from the one-loop quark-gluon vertex shown in Fig.lQd, following the usual one-loop 
PT procedure presented in section II. This part of the construction has been first carried out in |^9|] , where the 
process-independence of the IIq^ (g) was explicitly demonstrated. 

The second term, Tp^j.^^^ ((7,^1,^2), is Bose-symmetric with respect to pi <~> P2- As one can easily verify from Eq. 
(3.35), Tpl^ (g,pi,p2) is zero on shell, i.e. when contracting with the polarization tensors and using Pi = Pi = 0. 



Thus rp|^^j^(q,pi,p2) can be dropped when studying the one-loop process qq gg, exactly as rp^^^j^(q,pi,p2) was 

dropped in the tree-level case. As we we have seen in section IIIB, in the off-shell case, i.e. when Talviq,Pi,P2) is 
inserted into a two-loop quark-gluon vertex, the parts proportional to p^ and P2 pinch the internal quark propagator 
and give propagator-like contributions, whereas the parts proportional to pf and p| cancel exactly against analogous 
contributions from the rest of the graphs contributing to the two-loop quark-gluon vertex. Finally, Tp^^^^{q,pi,p2), 

is exactly the one-loop version of rp|^^j^((7,pi,p2): It is the one-loop three-point function involving one background 
(g) and two quantum (pi,P2) gluons as incoming fields, computed using the BFMFG Feynman rules. (Fig. 18). Notice 
that the special ghost structure characteristic of the BFM (Fig.lSe - Fig.lSh) emerges automatically, after following 
the procedure outlined above. 



(0) 



It is straightforward to show that rp^^^^(5,pi,p2) satisfies the following WI 

'?"r«„,.('7,pi.P2) = n«(pi) - n«(p2). 



(6.36) 



which is the exact one-loop analogue of the tree-level Ward identity of Eq (2.4); indeed the RHS is the difference of 
two conventional one-loop self-energies computed in the REG. 



In addition, we have for the action of pi^ on Fp^^ ((7,pi,p2) when pf 



pi 



0, is given by 



where Ila} (9) is given in Eq (^T4|) , and 



J3 



{k - pi)'T(,0) (P2, k, -k - P2) - (fc + P2).fca 



2B{q) + B{p^) 



(6.37) 



(6.38) 



P2 fcT^^Lpl?: fc + P2, -fc + Pi) - P2 ■ (fc - Pl)(2fc + P2 - Pl)a 
^ ■ B{pi)+B{p2) 

B{p)= I [dk]J,{p,k). 



with 



(6.39) 



(6.40) 



Eq. ( 6.37 ) is the one-loop analogue of Eq. ( 5.16 ). It is important to emphasize that rp)_^^j^((7,pi,p2) is not equal 

to the one-loop PT three-gluon vertex T'"al,u{q,Pi,P2) constructed in Q. Notice also that, unlike the tree-level case, 
now there will be a ^-channel ghost (we will not report its closed expression here) and that the t-channel has been 
modified also in order to achieve the PT one- loop rearrangement. 
Thus we have, 



14]F 



2s 



(16a) + (166) + (16d) -f (16e) 



(6.41) 



20 



with (we suppress a common factor g{'ya}f '^il) ) 



(i6&) = rW^.(g,Pi,P2), 

(16a) = nW(g)d(Q)r(?);jg,pi,p2), 
(16d) = 4°)^,(g,pi,p2)db2)n(i)'^(p2) , 

(16e) = r(?)^.(9,Pi,P2)d(pi)nW-(pi) (6.42) 



and so, together with the first equation in (5.17) 



rirrjr* = ((166) + (16d) + (16e) j . 7jr* + (16a) • tJ^* 



= |^(20a) + (206) + (20d) + (20/) + (20j) + (2O5) + (20/i) j + C2giuons | (4e) | (6.43) 

— (2) 

i.e. we recover the two-gluon Cutkosky cuts of ^apil^^Q ~ 1) (Fig-20), together with the corresponding two-gluon 

cuts of the IPR graph in Fig.4e, as we should. To account for the remaining two-ghost Cutkosky cuts in Fig. 20 and 
Fig.4e we need, in addition to Eq. ( |6.37| ), the following results 

pi . (16a) = -mW(g) +nW(g)d(g)(p2 -pi)>2,, 

Pi ■ (i6d) = i„,(g)n(i)^(p2) + inW(p2) , 

Pi-(16e) = 0. (6.44) 



Then we have 



Pi 



(16a) + (166) + (16d) + (16e) j = t^q) (^A^ + nW(p2)j 

+ (si'^+fli'^AQ)diq)iP2-Pir]p2. (6.45) 



From the above results one may immediately deduce the closed form of {Ai2^^)'^'' and ( S^}^^^'^^ ) appearing in Eq. 
(Hi: 



ab 



5[4]F{i2}y ^ gij-jrOdiq) +nW(g)d(g)(p2 -Pi)^j 

(Af'l^)f - ~^9m.r' (aW" + n(i)"(p2)^ (6.46) 

Then we have 

2^[4]F{i2}^[4]n2i} ^ 24^)(p2 - Pi), + 2mq)iP2 - PiUP2 - Plh 

= ((20z) + (20e) + (20m) + (20n) + (20g) + (20r)^ + C2ghosts | (4e) | (6.47) 
which concludes the proof of Eq ( |6.32D . 

B. Tree-level q{kx)q(k2) 3(^1)3(^2)3(^3), and the r*!,'^^^p((7,pi,p2,P3) 

Next we consider the tree-level amplitude for the process q{ki)q{k2) — > 5(^1)5(^2)5(^3), shown in Fig. 21 . This 
amplitude must be appropriately rearranged, and it will eventually furnish the imaginary parts corresponding to the 
three-particle Cutcosky cuts of &-'^^{q). 
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We start by presenting some general properties of this amplitude. It is straightforward to verify that T^l^^ satisfies 



the following identities: 



Pi - )p P^u + [^3 }u P3p , 



P'i'^'ilivp — {^3, 



p Pip + 1^3 )^ PZp ■ 
Ps-'S/Ji/p — Pa Iv Plp-^WS )p P2v 



(6. 



Bose symmetry imposes the following relations among the 5*3*"''^ amplitudes 

{sl''^)T^'^^{p,,p,,p,) = (5f 



(6.49) 



and 



P. 



be 



P 



(6.50) 



These sets of identities guarantee that all dependence on the unphysical four- vector Ua and the gauge parameter rj 
appearing in the polarization tensors will disappear from the final expression for ^3^' . 
The amplitude rearrangement is as follows (Fig. 21): 



irl!^^rp% - 9'{j'nd{q)Tf:;^;{q,Pi,P2,P3) 



(^j-^3]F\abc 



pup 

abc 
pvp 



(6.51) 



where 



rit;.p(9>Pl>P2,P3) = rT'd{p, +P2)T'ZpM^P3,Pl +P2)T%1''{P1,P2,-Pl -P2) 

+r'T''dih + fc3)r(,°L(9.P2,Pi +P3)r|5;(Pi,P3, -pi -ps) 

+ r-f^-dik, + h)TP^^Jq,p,,P2+P3)rl^^^P2,P3, -P2 -ps) 
I Y<(0)eabc 



(0) 



3.52) 



I t is e asy to verify now that r^^"'j'';(g,pi,p2,P3) is the analogue of r^^'^''j^(q,pi,p2) appearing in the first equation of 



(O)afcc/ 



(5.17) for the case of three on-shell gluons. In particular, we have that 



"4t;.'';(9:Pl:P2,P3) = -r^r'-'pldipi +P2)r<^Z(P^,P2,-pi - P2) 
-r''T''pld{h + k3)T^^^APl,P3,-pi -P3) 
_jea.jbc.p2 ^(^^ ^ fc3)r(0) (p2,P3, -P2 " P3) : 



(6.53) 



where we have used the well-known WI relating the bare three-and four-gluon vertices [Q. Eq. (6.53) is the analogue 
of Eq. (O). Clearly, when pl=pl=pl = 0, 



1 r>.' (Q,P1,P2,P3) =0, 



(6.54) 



exactly as happens in Eq. (2.4) 
In addition, 



Piq-[3]F^abc 
Pl\-^3s J pvp 

„M/'7-[3]\abc 
Pi \-^3t J pup 



(q[3\F {12} . abc ^ 
W3s )p P2v 



(q[3]F {13}. abc^ .,[3]_Fxq6c 
W3s lu P3p + (A3 . 



,'c.[3]{12}^abc„ , (q{3\{13}^abc^ fK[3]Fsabc 
[^3t )p P2u + {!^3t )u P3p - (,^3 



(6.55) 



where 



[Sf/^'^^)f^ = {^^nd{q) 



\mabc 



( r/ \mabc I / ry \mabc , { ry \mo 
[^l)ap + [^2)ap + (^3)ap 



(6.56) 
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with 



and 



{Z2) 



mabc 2 j^mbx £acx 

ap ^ 9 J J 



d{pi +P3)(2P2 +q)a{Pl +P3)p -9, 



T^-fp S{k2 + Ps) T^la + T^la S{ki + P3) -fpT^ 



U{k^)r'd{pi+P2)p"2 



Armed with the above relations, it is then straightforward to show that 



1/1 



2 V3! 



6 5. 



[3]F{12} _^ c.i:3]{12}^ ( ^[3]F{21} ^ ^[3]{21}y 



3s 



-yl'^^l +^['^1 +^[^1 



with 



3t J 



A 



[6] 



1 / 1 



2 V3! 



1 / 1 



^[3]Fj-l3]F* _ g^^[3]F{12} j ^^[3]F{21}j 



^ 2 V3! 
~ 2 I 3! 



T-[3]F^[3]* c. 
-'Ss ■'st ^0 '^38 



[3]F{12}\ /c.[3]{21}\ 



-[3] ^[3]* 
'3t -'3* 



,[3]F{12}\ /^[3]{21}y 



We are now in position to prove that, indeed. 



41^9) = ^3 nS(g) = C3 ni^j(g,eQ - 1) 



where is defined from exactly as in Eq. (mi) and Eq. (loil) 



To begin with, 



= (22a)c, + (22a),, + (226) + (22c) + (22.g) + (22£) , 



For the cuts involving ghosts we need t he cl osed form of S*. 



by virtue of the relations given in Eq. ( 6.4g| ). In particular, 



[3]F{21} 
3s 1 



which can be obtained from S". 



[3]F{12} 
3s 



of I 



{Z'Xpf^ = g^r^f'^^dip, +p,)p^,T%^^{q,p3,pi +P2) 



I \nabc 2 rnar rbcr 

(Z2)f3p ^9 J J 



d{p2 +P3) (2pi + q)f}iP2+P3)p - gpp 



{z'^TpT = g'r'T'^dip^ +P3) (2P2 + q)mp ■ 



Then it is straightforward to show that 

1 

2 
1 



ZiZ[ = (22/1), 

ZiZ^ = (22e)e, + (22A:), 
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1 , 

2^1^,3 


- (22e)ei 




2 1 


= (22/)c, 


+ (22j) , 


2^2 


= (22i),, 


+ (22.s) + (2271) + (22o) , 


\Z2Z'^ 


= (22<z) + 


(22m) , 


2 ^ 


= (22/)c. 




2^2 


= (22r) + 


(22p) , 


1 , 

2 3 3 


= (22«)e, 





(6.64) 

Thus, we have accounted for all three-particle Cutkosky cuts appearing in Fig. 22 . 

VII. DISCUSSION AND CONCLUSIONS 

In this paper we have presented the generalization of the PT to two-loops for the case of mass-less Yang-Mills 
theories. Two different, but complementary derivations have been presented. In the first derivation we have followed 
the diagrammatic approach employed in the one-loop case , and have shown how the PT properties are replicated 
in the next order. In the second we have pursued the dispersive PT construction established in and have shown 

that the resulting structures are consistent with those derived with the first method. We emphasize that throughout 
this entire analysis we have maintained a diagrammatic interpretation of the various contributions. In particular, no 
sub-integrations had to be carried out. This additional feature renders the method all the more powerful, because 
unitarity is manifest, and can be easily verified by means of the Cutkosky cuts. The combination of the two methods 
constrains significantly the PT construction presented here, and restricts severely any possible deviations from it. 
The reader should be able to recognize, for example, that any rearrangement of the (internal) vertices of two-loop 
box-diagrams (for example, of the so-called "H-diagram" discussed in the second paper of [0) cannot be reconciled 
with the arguments of the simultaneous two- and three-particle Cutkosky cuts presented in section VI. Whereas no 
strict proof has been given here that the PT construction developed in this paper is mathematically unique, we 
consider that as a very plausible possibility. Notice also that the appearance of characteristic one-loop structures 
inside the two-loop PT Green's functions suggests the onset of an iterative pattern, which may provide clues leading 
to the generalization of the PT algorithm to all orders in perturbation theory. 

The generalization of the two-loop PT construction to the case of Yang-Mills theories with spontaneous symmetry 
breaking (Higgs mechanism) in general, and the electroweak sector of the Standard Model in particular, should proceed 
precisely according to the methodology presented in this paper. Except for the additional book-keeping complications 
stemming from the presence of gauge-boson masses (modifications of WI, appearance of seagull and tadpole terms, 
diagrams with would-be Goldstone bosons), no additional conceptual obstacles are expected. 

The results of this paper clearly prove that the correspondence between PT and BFMFG persists at two- 
loops. Notice that this proof is based on an a-posteriori comparison with a result established through the systematic 
diagrammatic rearrangement of the physical S'-matrix computed in the renormalizable gauges, rather than on an a- 
priori formal derivation at the level of the BFM generating functional. It would be clearly important to reach a deeper 
understanding of what singles out the value = 1- One possibility would be to look for special properties of the 
BFM action at = 1 jH^ . In such a case one could choose to avoid the complications arising from renormalization, 
since the correspondence is valid also for the super-renormalizable 3-d QCD. 

Just as happened in the one-loop case, the two-loop PT self-energy defined here lends itself as an essential ingredient 
for the extension of the notion of the QCD effective charge 51 1^ to two-loops (for a thorough discussion of the 



one- loop case see |17[] ), since it has precisely the same properties as the corresponding QED quantity, i.e. the vacuum 
polarization of the photon. First of all, the two-loop PT self-energy captures the leading logarithms of the theory, i.e. 
the prefactor of t he lo garithm is the second coefficient of the QCD (3 function. Second, by virtue of the the QED like 
WI given in Eq. the combination aeff{q,fi) ~ g"^ (fi) A{q / fi) , where A(g/^) = [1 — U^^^q/fi) - n^^^fqV/i)]"^ is 

a renormalization group invariant quantity. Third, it has by construction the correct unitarity structure |14| . While 
the Oieff defined above appears as the obvious candidate, a detailed study needs be carried out in order to determine 
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whether or not there are any field-theoretical obstructions ||5^ which would prevent the realization of the two-loop 
construction corresponding to the QCD effective charge. 

In this context notice also that the construction presented here determines uniquely the constant term of the two- 
loop effective charge (within a given renormalization scheme) [ |5^ . To determine its actual value one should go beyond 
the two-loop calculation presented in and compute the Feynman diagrams of Fig. 9 keeping also constant terms 
. Knowledge of this constant term is important when the QCD effective charge is used in the field of renormalon 
calculus; in particular, it would be free of the ambiguities infesting the estimates of the renormalon contributions to 
physical observables |T7||5^ ]. 

In addition, it is well-known that if one was to compute the two-loop giuon self-energy in the context of the 
BFM keeping arbitrary, the resulting ^g-dependent term would be a constant, i.e. it would not affect the coefficient 
of the logarithm. While in such a case the residual gauge-dependence could be thought of as a renormalization-scheme 
ambiguity, i.e. it could be re-absorbed in the wave-function renormalization of the (background) gluon, this is not 
possible when the gauge fields are massive. In that case unitarity is even more constraining; as is known from the 
studies on the one- loop electro- weak effective charges |p"3|-|l5| , the gauge-dependence affects non-trivially the analytic 
structure of the answer, giving rise to unphysical thresholds. 

The calculations presented in section VI constitute the first dispersive derivation of the two-loop QCD /3 function. 
In this analysis we have made use of the one-to-one correspondence between the physical S'-matrix elements and the 
Cutkosky cuts of the two-loop PT self-energy. This construction involves a very particular combination of one-loop 
(section VIA) and tree-level graphs (section VIB); in addition to furnishing the correct P function coefficient, a subtle 
cancellation of infrared divergences also takes place: while both sets of graphs are infrared divergent, they combine 
to give a infrared finite answer. This can be directly inferred from the simple observation that the Cutkosky cutting 
procedure we have employed amounts finally to the determination of the imaginary part of a single logarithm, namely 
that of the two- loop self-energy; the latter is infrared finite [^,^. While the Cutkosky formalism furnishes an 
intuitive diagrammatic understanding and a valuable calculational short-cut, it would be interesting to reproduce 
the results of section IV without resorting to it. In particular one could study the precise cancellation mechanism 
of the infrared divergences using a proper infrared regularization scheme, and explicit expressions for the two- and 
three-gluon phase-space, which we have not needed here. In addition, it would be interesting to attempt a similar 
two- loop derivation using the formalism developed in and study possible connections. 

It has been often advocated that the non-perturbative QCD effects can be reliably captured at an inclusive level by 
means of an infrared finite quantity, which would constitute the extension of the perturbative QCD running coupling 
to low energy scales ]60| Early results by Cornwall based on the study of gauge invariant Schwinger-Dyson equations 
involving this quantity suggest that such a description can in fact be derived from first principles. According to 
this analysis, the self-interaction of gluons give rise to a dynamical gluon mass, while preserving at the same time the 
local gauge symmetry of the theory. The presence of the gluon mass saturates the running of the QCD coupling; so, 
instead of increasing indefinitely in the infrared as perturbation theory predicts, it "freezes" at a finite value P,|6ll. 
It would be interesting to revisit this issue in the light of the results derived in the present paper. For example, 
one could study the structure of the gauge-invariant Schwinger-Dyson equation for the PT gluon self-energy, and in 
particular the way the PT three-gluon vertex F^^^ enters in the PT gluon self-energy, using the two loop results as 
a guidance. In doing so one could hope to systematically improve on the analysis of where the gauge-technique 
ansatz for the vertex was used. In this context one may find it advantageous to rewrite the vertex F^'' appearing 
inside the two-loop PT gluon self-energy in terms of F^^^ ; one should then interpret the emerging residual terms as 
parts of the two-loop self-energy, even though they appear to be pinch-like i.e. once the PT self-energy has been 
fixed it may be recast into a different form, but no pieces should be re-assigned to vertices or boxes. 

Finally it would be interesting to pursue a connection with other field- or string-theoretical methods [|63[-|67t , either 
in order to acquire a more formal understanding of the PT, or in order to combine various attempts into a coherent 
framework. 
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Figure Captions 

Fig.l: Carrying out the fundamental vertex decomposition inside the non-abehan Feynman graph contributing to 
Fa ^ (a), gives rise to a genuine vertex (6) and a self-enery-Uke contribution (c). 

Fig. 2: The diagrammatic representation of PT one-loop self-energy H^^^, as the sum of the conventional self-energy 
^a/jj graphs (a) and (6), and the pinch contributions coming from the vertices (c). 

Fig.3: The PT one-loop vertex f^J^. 

Fig. 4: The one-particle reducible graphs before [(a), (6), (c), (d)] and after [(e), (/), (g), (h)] the PT rearrangement. 
G^^^ and G*^^) denote respectively the conventional and PT one-loop vertices with one-loop self-energy corrections to 
the external fermions included. 



Fig. 5: The PT rearrangement of typical one-particle reducible graph (a), giving rise to its PT counterpart (a^), 

P (c) and to Y^^ 



(2) (2) 

and to contributions to the first term of Fp (c) and to Yp (6). 



Fig. 6: The result of enforcing the PT decomposition on the external vertices of some of the two-loop Feynman 

(2) 

diagrams contributing the conventional two- loop quark-gluon vertex ' . 

Fig. 7: The result of enforcing the PT decomposition on the external vertices of some of the remaining two-loop 
vertex graphs. 

Fig. 8: The Feynman diagrams contributing to the conventional two-loop gluon self-energy H^^, in the gauges. 

— (2) 

Fig. 9: The Feynman diagrams contributing to the BFM two-loop gluon self-energy H^^. 

(2) 

Fig. 10: The one-loop counterterms contributing to the conventional two-loop gluon self-energy 11^^. 

Fig.ll: The one- loop counterterms necessary for the two-loop gluon self-energy H^^; they are identical to those 

~ (2) 

needed for the BFM two-loop gluon self-energy H^^. 



Fig. 12: The one-loop counterterms necessary to cancell the sub-divergences inside the conventional two-loop quark- 
luon vertex Fa ' . In the case of massive fer 
appropriate mass counterterm (not shown). 



gluon vertex Fa ' . In the case of massive fermions the wave- function counterterm K2^^ should be accompanied by the 



Fig. 13: The one-loop counterterms necessary to cancell the sub-divergences inside the PT two-loop quark-gluon 
vertex Fa ; they are identical to those neede for the BFM two- loop quark-gluon vertex F^ . 

Fig. 14: The fmidamcntal BRS-enforccd cancellation of s-channcl (a) and t-channel [{di) and (^2)] contributions, 
instrumental for the absorptive PT construction. Graph (6) gives rise to the correct ghost-structure. 

Fig.15: The Cutkosky cuts of the PT (and BFMFG) one-loop gluon self-energy. 

Fig. 16: The one-loop amplitude for the process q{ki)q{k2) g{pi)g{P2), after the PT rearrangement. 
Fig. 17: The Feynman diagrans contributing to the conventional one-loop three-gluon vertex Ta^u- 
Fig. 18: The diagrammatic representation of the one-loop three-gluon vertex F^^^^^. 
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Fig. 19: Some of the one-loop t-channel graphs contributing to q{k\)q{k2) — > g{pi)g{p2)- 

Fig. 20: The two-particle Cutkosky cuts of the PT (and BFMFG) two- loop gluon self-energy. We have used the 
same labelling of individual diagrams as in Fig. 9. The two upper (lower) rows show graphs where two gluon (ghost) 
lines have been cut. 

Fig. 21: The tree-level graphs contributing to the process q{ki)q{k2) — > g{pi)g{p2)g{Pi), after the PT rearrange- 
ment. 

Fig. 22: The three-particle Cutkosky cuts of the PT (and BFMFG) two-loop gluon self-energy. We have used the 
same labelling of individual diagrams as in Fig. 9. The first fivc^ graphs have three-gluon cuts, the next two have 
two-gluon-one-ghost cuts, while the remaining ones have one-gluon-two-ghost cuts. 
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